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L'invention est relative k un proc^dd pour extraire un horizon g€oIogique et 
des propridtds associees de leur image ddduite de donn^es sismiques. 

L'invention est dgalement relative k un dispositif pour extraire un horizon 
geologique et des propridt^s associ^ de leur image d€duite de donndes sismiques. 

L'invention est enfin relative k un produit-programme d'ordinateur pour 
extraire un horizon g€ologique et des propri^t€s associees de leur image d^duite de 
donnas sismiques. 

On connalt des procdd^s d'extraction d'horizon et de leurs aflributs sismiques. 
Dans ces procddfe connus, on construit une matrice d'attribut sismique 
tridimensionnelle de m8me taille et de m6me fchantiUonnage que la matrice 
sismique tridimensionnelle inidale. 

Pour attribuer k chaque point de rhorizon une valeur d'un attribut sismique, 
on pent proc^der soit par interpolation des noeuds environnants de la matrice 
d'attribut sismique tridimensionnelle ou choisir la valeur stock^e dans le noeud le plus,- 
proche de la matrice d'attribut sismique tridimensionnelle. 

En cas d'interpolation, le pr6-traitement correspondant n€cessite un temps de ; 
traitement important et une quantity de memoiie disponible trfes importante ; c'est 
pourquoi il est gdn^ralement n&essaire d'effecmer ce pre-traitement sur un 
ordinateur puissant. 

Dans le cas du choix de I'attribut du nceud le plus proche, ce pr6-traitement 
n'est pas necessaire, et la determination de I'attribut ne pent Stre effectude qu'^ des 
points voisins des points du rdseau de la matrice sismique tridimensionnelle 
origiuale. 

Un inconvenient majeur de I'art anterieur est d'introduire des ddfauts ou 
artefacts qui induisent des eireurs verticales sur les horizons et les attributs sismiques 
et conduisent k une mauvaise interpretation g&logique des mesures sismiques. 

Un premier but de l'invention est d'amdliorer la technique connue en 
proposant un calcul direct et plus exact, en minimisant les risques d'erreurs possibles 
sur les attributs en tons points de I'horizon. 
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Un deuxifeme but de rinvention est de supprimer les erreurs verticales 
susceptibles de conduire a une mauvaise interpretation geologlque. 

L*invention a pour objet un procedd pour extraire un horizon g^ologique et 
des propri6t6s assocides d'une representation sismique, dans lequel on construit une 
5 fonction continue Sij,k(t) par interpolation ou approximation des traces sisniiques 
discretes d'une matrice sisnndque multidimensionnelle, ladite fonction etant ddsignde 
coDcune "trace sismique locale continue", comportant les Stapes suivantes : 

a) ^ utiliser comma decalage optimal de deux traces sismiques locales 
continues voisines, la valeur de decalage rendant maximale leur fonction de 

10 correlation ; 

b) . retenir coimne voisinage conditionnel d'une trace sismique locale continue 
"centrale" Sij,k(t) le sous-voisinage consistant en des traces voisines Spq,k(t) 
correspondant k des ddcalages optimaux associes k des correlations Rij,pq,k(h) 
superieures k un seuil predetermine compris entre 0 et 1 ; 

15 c). construire une matrice bidimensionnelle d'extraction destinee k Stre 

remplie par des valeurs extraites ; 

d) . choisir un point «graine» P(ij,t) et determiner le point P(i,j,k) 
verticalement le plus proche ; 

e) . estimer les proprietes associees du voisinage conditionnel et remplir la 
20 matrice bidimensionnelle d'extraction par des proprietes decalees par translation de la 

variable courante (t) de la valeur du decalage optimal (hij,pq,k) correspondant au point 
P(i j,k) verticalement le plus proche. 

Selon d'autres caracteristiques de Tinvention : 

- on remplace successiveraent le contenu de la matrice bidimensionnelle de 
25 retape c) au cours des iterations successives d'extraction, 

- k chaque iteration, on utilise comme nouveaux points «graines» de retape 
d) tons les points correspondent k des emplacements de la matrice bidimensionnelle 
de retape c) remplis k Titeration precedente, 

- la propriete de subsurface extraite est Tamplitude sismique refiechie et 
30 captee par des geophones. 
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- on visualise sur un ^cran de visualisation des surfaces horizon peintes avec 
des couleurs correspondant a un codage des propri^tds extraites. 

L'invention est ^galement relative h un dispositif pour la raise en oeuvre du 
procddd selon l'invention, comportant des moyens pour utiliser comme ddcalage 

5 optimal de deux braces sismiques locales continues voisines, la valeur de d&alage 
rendant maximale leur fonction de correlation, des moyens pour retenir comme 
voisinage conditionnel d'une trace sismique locale continue "centrale" Sij.ic(t) le sous- 
voisinage consistant en des traces voisines Spq,k(t) con-espondant k des d&alages 
optimaux associes k des correlations Ry,pq.k(h) supdrieures h. un seuil prddaermine 

10 compris entre 0 et 1, des moyens pour construire une matrice bidimensionnelle 
d'extraction destinde k Stre remplie par des valeurs extraites, des moyens pour choisir 
un point «graine» P(i,j.t) et determiner le point P(io,k) verticalement le plus proche 
et des moyens pour estimer ies propri6tes associees du voisinage conditionnel et 
remplir la matrice bidimensionnelle d'extraction pai- des proprietes d^caiees par 

15 translation de la variable courante (t) de la valeur du d^calage optimal conespondant 
au point P(i j,k) verticalement le plus proche. 

Selon d'autres caracteristiques de l'invention : 

- le dispositif comporte des moyens de memorisation et des moyens de 
visualisation de paramfetres sismiques determines k I'aide du precede selon 

20 l'invention. 

L'invention a enfin pour objet un produit-programme d'ordinateur, 
comportant des elements de code de programme pour executer Ies etapes d'un 
procede selon l'invention, loi-sque ledit programme est execute par un ordinateur. 

L'invention sera mieux comprise grSce k la description qui va suivre donnee k 
25 titre d'exemple non limitatif en reference aux dessins annexes dans lesquels : 

- La figure 1 represente schematiquement une matrice sismique 
tridimensionnelle et une tr^e sismique locale continue. 

- La figure 2 represente schematiquement, un exemple de voisinage local Nij,k 
consistant en un ensemble de trace sismique locale continues Spqjt(t) voisines de la 

30 h-ace centrale de reference Sijjc(t) elle mSme situee a la verticale du geophone Gij et 
interpolant les amplitudes sismiques autour de t=k. 
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- La figure 3 repr&ente sch^matiquement, une section verticale de voisinage 
local conditionnel Ny^kCr) d'une trace sismique locale continue de r^f&ence Sij,k(t). 

- La figure 4 repr^gsente sc^matiquem^l^junejection verticale d^ine matrice 

sismique tridimensionnelle avec un d^calage optimal hij,pq,k et un voisinage 
5 conditionnel Nij,k(r).. 

- La figure 5 repi'dsente schdmatiquement, une section verticale d'une matrice 
sismique tridimensionnelle analogue h. la figure 4 avec presence d*une faille ou 
discontinuity 

- La figure 6 repr^ente scMmatiquement un organigranmie fonctionnel d'un 
10 proced€ selon Tinvention. 

- La figure 7 reprdsente schematiquement, une vue partielle agrandie d'une 
noatrice sismique tridimensionnelle pr^sentant une discontinuite et illustrant la mise 
en oeuvre du proc€d6 selon Tinvention. 

En r^f^rence aux figures 1 ?t 5, les ^dments identiques ou fonctionneilement 
15 Equivalents sont d^ign^s ou rcp6x6s de mani^re identique. 

Sur la figure 1 , une matrice sismique tridimensionnelle est obtenue par releve 
de mesures enregistr^es par des g€ophones Gy disposes sur un r^seau x,y en des 
points de coordonn^es i, j. Le relev6 des mesures dchantiUonndes dans le temps est 
represents suivant un axe t descendant reprdsentatif de la profondeur ou d'une 
20 verticale descendante k partir de la surface de la terre ou de la mer. Les mesures sont 
caractSiisSes par leur amplitude, par exemple une amplitude relevSe par le gdophone 
Gy au temps ou h la profondeur d'&hantillonnage tk. La mesure discrete effectuSe par 
le gSophone Gij au temps ou k la profondeur tk est appel^ amplitude sismique Syje- 

L'ensemble des amplitudes sismiques correspondant k un geophone Gy de 
25 coordonnees i, j est une matrice unidimensionnelle (Syi, Sya,...., Syk,..-, S^n) 
appelSes trace sismique discrete, car cette matrice unidimensionnelle correspond a la 
tiace selon le point de coordonnees horizontales i, j de la matrice sismique 
tridimensionnelle obtenue par mesures sismiques. 

Uaxe vertical t orient^ selon une verticale descendante dSsigne usuellement le 
30 temps, mais pent dgalement Stre traits pour repr&enter une profondeur k partir de la 
surface. 
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L'invention conceme aussi bien rapplication k une troisifeme coordonnde t 
representative du temps, qu*^ une troisifeme coordonnde t representative de la 
profondeur. 

A partir de la trace sismique discrete situ6e k la verticale d'un g^Dphone Gij 
5 on definit, par interpolation ou approximation des valeurs discretes autour de t=tie=k, 
une fonction continue SyjcCt) qui est d&ignde comme "trace sismique locale 
continue", Les m6thodes d'approximation ou d'interpolation de valexirs discrfetes pour 
engendrer une fonction continue sont nombreuses, et comprennent notamment les 
interpolations ou approximations polynomiales, ainsi que les interpolations ou 
10 approximations trigonomdtriques polynomiales. 

Toute autre variante d*interpolation ou d'approximation foumissant une 
fonction continue pent 6galement 6tre appliquee k la pr&ente invention pour fourair 
une "trace sismique locale continue". 

Sur la figure 2, plusieurs traces sismiques locales continues ddfmissent un 
15 voisinage d'une trace sismique locale continue "centrale" Sij,k(t) de reference. Le 
voisinage d'une trace sismique locale continue Sy^kCt) est defini comme Tensemble 
des traces sismiques locales continues dont les indices spatiaux horizontaxix p,q sont 
voisins des indices spatiaux horizontaux i, j de la trace sismique locale continue de 
reference. 

20 A titre d'exemple, les coordonnees spatiales horizontales p,q correspondant au 

geophone Gpq sont voisines des coordonndes spatiales horizontales i, j correspondant 
au geophone Gij si les valeurs absolues des differences i-p et j-q sont inferieures k des 
entiers donnas, par exemple k 2. 

Dans ce cas, comme represents kla figure 2, la trace sisnuque locale continue 
25 Sij,k(t) est associee k huit traces sismiques locales continues voisines entourant la 
trace sismique locale continue "centrale" Sij.k(t). 

Dans le cas de traces sismiques locales continues produites par des mesures 
sismiques, le profil des horizons geologiques introduit des decalages verticaux entre 
des traces sismiques locales continues voisines. Dans le but de determiner des 
30 relations entre deux traces sismiques continues locales voisines Sij,K(t) et Spq,k(t) 
centi-ees sur la m6me coordonnee verticale d'echantillonnage t=tk=k et correspondant 
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h, des coordonnees spatiales i, j et p,q differentes, on calcule la fonction de correlation 
Rij,pq.k(h) des deux traces sismiques locales continues. 

La fonction de correlation de deux traces sisnaiques locales continues 
voisines est obtenue par la formule suivante 

5 R Ml) ^ ^^^''^'^ 

VCmK(0).C^.^..(0) 

oh le numerateur correspond k la fonction de covariance de S5j,k(t) et Spqjc(t) 

tk+A 

obtenue par I'expression suivante Cij,pqjc(h) = JSy (t) , Sj^ „ (t + h) . dt . 

Dans cette int^grale d^finissant Cij,pq.k(h) , le paramfetre A d^finit une «fen8tre 
10 verticale d'investigation » autour de t?=tic=k. 

Par exemple, si les traces sismiques locales continues Si,\k(t) et Spq,k(t) sent 
des polyn6mes trigonometriques de la forme suivante interpolant les donnas 
sismiques 

m 

Sij.k(t) = ^s^^'-c^sCsO) t) -h bf ^'.sinCsfi) t) 

5=1 

m 

15 Sij^(t) = ^» * -cosCsfi) t) + b*'' .sin(sa) t) 



alors on peut montter que la fonction de covariance Qj,pqjc(h) est elle mSme un 
polyn6me trigonom^trique de la fonne smvante oii les coefficients As''''^* et Bs^'^^ 
d^pendait des coefficients as*"*, 2is^, bj^ et bs^* : 



Cij.M.k(h) = XAiJ-'^\cos(sfi)h) + Bii-W'\sin(sfi)h) 

Un r&ultat mathdmatique classique et connu est que la fonction de correlation 
R««uc(h) = ^u>M>»^ traduit une similarity des fonctions correlees lorsque 

25 cette fonction de corr^ation se rapproche de 1 . 



L'etude des fonctions de corr^ation des traces sismiques locales continues 
Sij,k(t) et Spq^(t) permet de d^finir conraie decalage optimal la valeur hy.pq^k de h 
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correspondant au maximum de la fonction de correlation Rij,pq,k(h), c'est-k-dire au 
maximum de la fonction de corr^ation le plus proche de 1. Contrairement a Tart 
existant has6 sur une formulation discrete de la fonction de correlation Rij.pq,k(h), 
I'utilisation d'une formulation continue des traces sismiques et done de la fonction 
5 de correlation permet d'obtenir un ddcalage optimum qui n'est pas contraint a etre un 
multiple entier du pas d'echantiUonnage suivant Taxe vertical correspondant h, la 
variable L 

La definition du ddcalage optimal hij,pq,k de deux traces sismiques locales 
continues voisines permet d'obtenir une premiere approximation des horizons passant 
10 par cette trace sismique locale continue de reference Sij,k(t), comme represente aux 
figures 3^5. Contrairement a I'art existant, le fait d'utiliser des ddcalages optimaux 
qui ne sont pas contraint k etre un multiple entier du pas d'echantillonnage suivant 
Taxe vertical perniet d'eviter des erreurs connues sous le nom de «aliasing». 

Sur ces figures 3 a 5, on d^finit un voisinage conditionnel Nij.k(r) de la trace 
15 sismique continue locale de reference Sij,k(t) qui est un sous-voisinage du voisinage 
initial des traces sisnoiques locales continues retenues pour le calcul de coixdlation et 
de definition des ddcalages optimaux. Le voisinage conditionnel Nij,k(r) est choisi de 
manifere que, pour toute trace sismique continue locale Spq,k(t) appartenant a Nij,k(r), 
le ddcalage optimal hij,pq,k correspondant est tel que la correlation Rij,pq.k(hij.pqjc) entre 
20 Sij,k(t) et Spq,k(t) est superieure k un seuil predetermine r compris enti-e 0 et 1 . 

Comme on le voit sur la figure 5, la fixation du seuil r permet egalement 
d'englober dans un voisinage conditionnel une discontinuite ou une faille geologique, 
ce qui constitue un avantage important par rapport k Tart anterieur. 

Les dispositions precitees fournissent ainsi une modeiisation continue 
25 permettant la mise en oeuvre de Tinvention. 

En reference k la figure 6, un precede selon Finvention comporte une etape de 
modeiisation nuraerique 100 pour definir des traces sismiques locales continues, des 
decalages optimaux et des voisinages conditionnels, conmie decrit en reference aux 
figures 1 5. 

3^ A retape 101, on definit une mabice bidimensionnelle dont les indices de 

lignes et de colonnes correspondent aux coordonnees des geophones de mesures 
sismiques. 




A Tetape 102, on choisit un point «graine» P(ij,t) ou premier point 
d'extiaction. 

A r^tage 103, on detern xine le Jioe^d de la matrice sismique 

tridimensionnelle le plus proche du point «graine» P(i j,t). 

5 Ce noeud a pour coordonn^es spatiales i j et pour coordonnde temporelle ou 

de profondeur la coordonnfe tk=k la plus voisine de la coordonn^e t du point 
«graine» P(i j>t) choisi Jt T^tape 102, 

A rdtape 104, on s^lectionne tout d*abord les traces sismiques continues 
locales Spqjc(t) appartenant au voisinage conditionnel Ny,k(r) defini k I'^tape 100 et 
10 telles que les indices (p,q) correspondent a un emplacement vide de la matrice 
bidimentionnelle definie a T^tape 101. 

Ensuite, pour chaque trace sismique continue locale Spq,ic(t) ainsi selectionnee, 
on attribue au point P(p,q,t) la valeur P(p,q,t + hij,pqjc) , oft hij,pq,k est le decalage 
optimal du point P(p,q,t) par rapport au point voisin PCiJJk:) determine & Tetape 103. 
15 Ce point P(p,q,t + hy^pq^c) est consid^r^ comme un nouveau point situe sur Thorizon 
et est stocks dans Templacement (p,q) de la matrice bidimensionnelle d6fini a Tfitape 
101. 

Le precede qui vient d'etre d€crit est destine k determiner un ensemble de 
points appartenant a Thorizon passant par un point P(ij,t) , mais est egalement 
20 applicable pour attacher ^ ces points, au fur et k mesure de leur determination, un 
ensemble de propri^tes caracterisant la nature physique des terrains au voisinage de 
rhorizon ainsi determine. 

A cet effet, on peut prevoir d'attacher k chaque point P(p,q,t + hy^pq^k) de la 
matrice bidimensionnelle definie h Tetape 101 un ensemble de proprietes appeld 
25 "attributs sismiques", 

Ces propridtds sont gen^ralement representees pai* des coulem-s differentes sur 
les visualisations d'ecran, de maniere k permettre une visualisation rapide et sans 
risque d'erreur des proprietes associees kun horizon geologique. 

En pratique, on effectue des etapes sinoilaires aux etapes decrites en reference 
30 aux etapes 100 k 104 et Von construit les attributs du point P(p,q,t + hy^pq^k) k Taide 
de requation de la trace sismique continue locale Spq,k(t) ; par exemple on construit 
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ainsi un attribut enveloppe. un attribut phase ou un attribut frequence, de maniere 
connue en soi, par exemple d'aprfes I'enseignement du document de TANER M.T., 
KOEHLER F..SHER1FF R.E., (1979), Complex seismic trace analysis. GEOPHYSICS.' 
volume 44, n° 6, pages 1041 h. 1063. 

5 Cette mdthode permet ainsi d'afficher I'attribut sismique sur l-hcrizon extrait 

en le "peignant" avec les couleurs conespondantes. 

Le proc^dd venant d'etre d6crit permet ainsi de definir des vateurs d'extraction 
d'un horizon pr^daermine correspondant k une donne'e r^elle, mais uniquement dans 
un voisinage conditionnel du point P(i,jjc) le plus proche du point «graine» P(io,t) 
10 choisi. 

Pour propager le procdd^ d'extraction h I'ensemble de la matrice 
bidimensionnelle d^finle k V6tape 101, on construit tout d'abord un ensemble de 
nouvelles graines constitu^ par les nouveaux points P(p,q,t + hij,pq.k) situe sur 
rhorizon qui viennent d'etre d€termin€s et est stocks dans les emplacements (p.q) 
correspondants de la matrice bidimensionnelle defmi k I'e'tape 101. Chacun de ses 
pomts est memorise k I'^tape 106. pour etre r^injectd k I'^tape 102, en y servant de 
nouveau point «graine» . 

Ceci propage le traitement successivement k I'ensemble de la matrice 
bldimensionneUe correspondant ame gdophones. 

Lorsque qu'aucune nouvelle graine ne pent 6tre trouv^e, le proc€d^ est 
orients k I'etape 106 vers une dtape 107 de visualisation sur dcran des horizons ou 
des attributs ou propri6t& assocides d'horizon geologique, de manifere k pennettre 
une visualisation en couleur ou en niveaux de gris des propridtds assocides k 
I'horizon geologique extrait correspondant au temps i^l t. 

25 Ainsi, ^ la fin de I'algorithme d^fini en r6f6reac& k la figure 6, la matrice 

bidimensionnelle de I'^tape 101. contient une plurality de points correspondant au 
meme horizon que le point initial «graine» choisit P(iJ,t). 

L'ensemble des points pent etre i^ordonnd avec ses voisins, confonn^ment k 
la figure 7 pour definir un maiUage recouvrant entiferement I'horizon contenant le 
30 point «graine» P(ij.t). ce maillage permet ainsi de repr^enter I'horizon extrait 
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corame une surface localement continue sur I'&ran de visualisation d*une station de 
travail. 

_Le procede selon rinvention est de pr6f€rence mis en oeuvre sur un dispositif 
comportant des moyens appropri^s a la mise en o&uvre des Stapes decrites en 
5 r6f6rence k la jHgure 6, 

En particulier, un dispositif selon Tinvention comporte des moyens de 
m^naorisation necessaires aux calculs successifs et les naoyens de visualisation 
necessaires & T^tape 107. 

De pr^f^rence, un dispositif selon Tinvention est un dispositif programmable 
10 commande par un produit-progranome d'ordinateur comportant des 6l6mtnts de code 
de programme pour ex€cuter les Stapes du iproc6d6 d6cnt en r6fdrence la figure 6. 

L'invention decrite en reference k plusieurs objets particuliers n'y est 
nullement limitee, mais couvre au contraire toute modification de forme et toute 
vaiiante de realisation dans le cadre et Tesprit de I'invention. 
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REVBNDICATIONS 

1. Procede pour extraire un horizon g^ologique et des propri6t6s associees 
d*une representation sismique, dans lequel on construit une fonction continue Sij,k(t) 
par interpolation ou approximation des traces sismiques discretes d'une matrice 

5 sismique multidimensionnelle, ladite fonction 6tant d€sign€e comme "trace sisniique 
locale continue", comportant les Stapes suivantes : 

a) , utiliser comme decalage (vertical) optimal de deux traces sismiques 
locales continues voisines, la valeur de decalage rendant maximale leur fonction de 
conelation, ce decalage optimal n'etant pas contraint k gtre un multiple entier du pas 

10 d'echantillonnage vertical ; 

b) retenir comme voisinage conditionnel d'une trace sismique locale continue 
"centrale" Sij.k(t) le sous-voisinage consistant en des traces voisines Spq.it(t) 
correspondant & des ddcalages optimaux hj^pq^ associ& h des correlations 
Rij.pq.k(hij.pq,k) superieuics k un seuil pr€d6termind compris entre 0 et 1 ; 

15 c). construire une matrice bidimensionnelle d'extraction destinee a etre 

remplie par des points extraits appartenant au mSme horizon que celui passant par le 
point graine ; 

d). choisir un point «graine» P(i,j,t) et determiner le point P(i j,k). de la 
matrice sisnndque tridimensionnelle verticalement le plus proche ; 

20 e). estimer les proprietes associees du voisinage conditionnel et remplir la 

matrice bidimensionnelle d'extraction par des proprietes decalees par translation de la 
variable courante (t) de la valeur du decalage optimal (hy^pqjc) correspondant au point 
(i, j, k) verticalement le plus proche. 

2. Froc6d6 selon la revendication 1, dans lequel on utilise comme point 
25 «graine» de Tdtape d) tous les points nouvellement stockes dans la matrice 

bidimensionnelle de Tetape b) et non encore utilises comume point graine. 

3- Precede selon la revendication 1 ou 2, dans lequel on remplit 
successivement le contenu de la matrice bidimensioimelle de Tdtape c) au cours des 
iterations successives d'extraction. 




4. Proc6d6 selon la revendication 1 ou 2, dans lequel on remplace 
successivement le contenu de la matrice bidimensionnelle de Tdtape c) par ime 
moyenne des contenus successifs au cours des iterations successives d'extraction. 

5. Froc6d6 selon la iievendication 1 ou 2, dans lequel la propriete de 
5 subsurface extraite est un ensemble d'attributs sismiques calcules en chaque point 

extiait sur Thorizon passant par le point graine, le calcul de ces attributs 6tant lui 
naSme r^alis^ simultan6ment a T extraction de ces points. 

6. Procdde selon la revendication 1 ou la revendication 5, dans lequel on 
visualise sur un €cran de visualisation les attributs sismiques « peints » sur Thorizon 

10 extrait. 

7. Dispositif pour la mise en oeuvre du procede selon Tune quelconque des 
revendications 1^4, comportant des moyens pour utiliser conune decalage optimal 
de deux traces sismiques locales continues voisines, la valeur de ddcalage rendant 
maximale leur fonction de correlation, des moyens pour retenir comme voisinage 

15 conditionnel d'une trace sismique locale continue "centrale" de reference Sij,k(t) le 
sous-voisinage consistant en des traces voisines SpqjcCt) correspondant h des 
d^calages optimaux hij,pq,k associes k des correlations Rij,pq.k(hij,pq^) sup6rieures a un 
seuil predetermine compris entre 0 et 1, des moyens pour construire une matrice 
bidimensionnelle d'extraction destine k etre remplie par des valeurs extraites, des 

20 moyens pour choisir un point «graine» P(ij,t) et determiner le point P(ij,k) 
verticalement le plus proche et des moyens pour estimer les proprietes assocides du 
voisinage conditionnel et remplir la matrice bidimensionnelle d'extraction par des 
propri6t6s decal€es par translation de la variable courante (t) de la valeur du decalage 
optimal (hij,pq^) correspondant au point (i, j, k) verticalement le plus proche. 

25 8» Dispositif selon la revendication 7, comportant des moyens de 

memorisation et des moyens de visualisation de paramfetres sismiques determines k 
Taide du procede selon Ihine quelconque des revendications 1 & 6. 

9. Produit-progranome d'ordinateur, comiportant des elements de code de 
programme pour executer les etapes du procede selon Tune quelconque des 
30 revendications 1^5, lorsque ledit programme est execute par xm ordinateur. 
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10. Produit-programme d'ordinateur, comportant des 616ments de code de 
programme pour executer les etapes du procede selon la revendication 6, lorsque 
ledit programme est execute par un ordinateur. 




FIG. 2 




FIG. 5 
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